Tunneling of polarized fermions in 3D double wells 
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We study the tunneling of a spin polarized Fermi gas in a three-dimensional double well potential, 
focusing on the time dynamics starting from an initial state in which there is an imbalance in the 
number of particles in the two wells. Although fermions in different doublets of the double well 
tunnel with different frequencies, we point out that (incoherent) oscillations of a large number of 
particles can arise, as a consequence of the presence of transverse degrees of freedom. Estimates 
of the doublet structure and of the occupation of transverse eigenstates for a realistic experimental 
setup are provided. 
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I. INTRODUCTION 



Tunneling of particles through wells and barriers is a distinctive property of quantum mechanics, commonly exploited 
in the realization of solid-state devices [lj, |2| . The experimental realization and manipulation of ultracold atoms [3|, |J] 
in double- and many- well potentials makes possible to study tunneling dynamics of fermions and bosons in an highly 
controllable setup, in which it is possible to tune the geometrical properties of the wells, and then the tunneling rate 

£rv, of the particles. For ultracold bosons, the coherent Josephson oscillations of a Bose-Einstein condensate in a double 
well potential were discussed [a, Q and subsequently experimentally observed [3, Q ; the dynamics of bosonic squeezed 
states has been as well investigated [9J]. For large enough barriers between the two wells, a two-mode ansatz can 

d ' be used to describe the tunneling dynamics: for bosonic condensates, the two-mode equations can be mapped in a 

3 ^ non-rigid pendulum [5J and can also generalized to many- well potentials [10 . Ill| . 

On the other hand, the study of ultracold Fermi gases has known a great advance through the last years |12Hl4| and 
the impressive experimental progresses in their control make Fermi gases very promising to realize ultracold fermionic 
junctions and to study tunneling phenomena. Due to the possibility to use optical lattices |15l - [l8j and to tune the 
interspecies interaction through Feshbach resonances [3J, lj] , the fermionic tunneling can be studied in situations having 
a direct counterpart in solid-state tunnel junctions. When the Fermi gases in different wells are superfluid (due 
to attractive interaction among species) one has the ultracold fermionic equivalent of a superconductor-insulator- 
superconductor junction [lj], having coherent tunneling: theoretical studies of tunneling of fermionic superfluids 

O • through barriers [19-21] and in double- and many- well potentials [22 H28[ have been reported in literature, as well as 
the study of the internal Josephson tunneling between different species [29| (see more references in |21|). 

— , . At variance, one can also realize the ultracold atomic counterpart of one or more normal tunnel junctions [301 ] 

^ ' when the temperature of the Fermi mixture in a double- or many- well potential is above the superfluid critical 
temperature (in presence of attractive interactions) or, even at zero temperature, when the Fermi gases in the wells 
are in the normal state. The latter situation can be obtained polarizing above a critical threshold a two-component 
Fermi mixture [3l|, |32| or simply using spin polarized fermions [33J, |3J]. In [331 ] the center-of-mass motion of a 
polarized Fermi gas in a combined periodic and harmonic potential was theoretically and experimentally investigated, 
showing an insulating regime when the Fermi energy lies into the bandgap of the lattice: working in tight-binding 
approximation, the dynamics of the polarized gas was studied solving the classical Liouville equation, allowing for the 
characterization of the different regimes of the center-of-mass dynamics of the 3D Fermi gas |33|, |35J] . 

In this paper we study the tunneling dynamics of an ideal (spin polarized) fermionic gas in double well potentials. 
Of course, unlike condensed Bose gases where a large number of particles are in the same state moving coherently, for 
polarized fermions the dynamics displays in general an incoherent motion of the particles in the external potential. The 
situation is analogous to the incoherent tunneling of electrons between normal metals in the approximation where they 
are considered free particles: the application of a constant voltage creates a difference in the chemical potential and 
the more energetic electrons can tunnel through the barrier. For Fermi gases an imbalance in the chemical potential 
can be either obtained by making asymmetric the double well potential (with the two wells having different potential 
minima) or creating an imbalance in the number of particles between the two wells of the system at the initial time: 
in this paper we will focus on the latter situation. Because of Pauli principle, fermions occupy all the doublets of 
the double well potential up to the Fermi energy: since in general fermions in different doublets tunnel with different 
frequencies, then a dephasing in the current flowing among the wells may arise. This is what happens in the ID case 
[34J : the transverse degrees of freedom are frozen, since the distances between different doublets of the ID double 
well potential are much smaller than the confining frequencies in the transverse directions. The tunneling dynamics 
then shows strongly aperiodic spatio-temporal patterns [34J : for large number of fermions, current oscillations are 
practically washed out. Rabi oscillations of a degenerate fermionic gas in a double well potential has been studied in 
[36| in relation to the possibility to exploit them for the interferometric measurement of external forces at micrometer 



length scales. 

However, for a relatively large number of fermions (say N > 10 3 ), the validity of ID limit requires very large confining 
transverse frequencies oj±: indeed, as we will discuss in Section II, if one has an energy barrier Vq between th e wells 
and the minima of the wells are at distance A/2, then the energy difference between the doublets is ~ UkyVoJm 
(where k — 2ir/X and m the mass of the fermionic atom). The condition for which the transverse degrees of freedom 
are frozen is then Nhky/Vo/m < fruj_. For a realistic double well potential [7| one has Vo/h ~ 1kHz and A ~ 2/j,m: 
with m the mass of potassium ( 40 _ftT) atoms, one gets lu±/(2ttN) > 1kHz, corresponding to transverse frequencies 
u)±/2ir > 1000kHz. One then sees that for typical transverse frequencies one has to take into account the transverse 
degrees of freedom. 

In the following we consider a realistic 3D double well, obtained superimposing an harmonic confinement with a 
periodic potential: properly tuning the parameters, one can have two wells much more populated than the others [7|. 
We provide estimates of the tunneling rates and of the occupation of the transverse degrees of freedom. In particular, 
we s how t hat it is realistically possible to have only the first doublet occupied: this corresponds to the condition 
hky^Vo/m > Nhio±, showing the role played by the transverse degrees of freedom in providing a reservoir to store 
particles and resulting in oscillations of a large number of particles with the same frequency. We stress that these 
are incoherent single-particle oscillations. Then, if the ID condition is violated the simple observation of a sinusoidal 
current cannot discriminate between incoherent and coherent tunneling dynamics (the latter being expected when the 
weakly coupled Fermi gases are superfiuids). 

The plan of the paper is as follows: in Section II we analyze the structure of levels for a 3D Fermi gas in a double 
well potential obtained by the superposition of an harmonic confinement and an ID optical lattice. In this Section we 
also provide estimates for the number of doublets and the occupation of transverse eigenstates in a realistic setup. In 
Section III we discuss the equations of motion for the dynamics in a 3D double well, and we show how the ID limit, 
studied in [34[ , is retrieved. In Section IV we study the dynamics of the Fermi gas in the 3D double well potential 
described in Section II, showing how oscillations characterized by a single frequency can arise and studying deviations 
due to the occupation of higher doublets. Finally in Section V we present our conclusions. 

II. POLARIZED FERMIONS IN A 3D DOUBLE WELL POTENTIAL 

In this Section we study the structure of levels and their filling by a Fermi gas in a 3D double well potential. To 
fix the ideas, we focus on the double well potential obtained by superimposing an harmonic confinement and a ID 
optical lattice: when the energy of two minima of the periodic potential is significantly lower than the energies of 
other minima, due to the presence of the harmonic trap, then these wells are much more populated than the others 
and one practically has a double well. This way to realize the double well has been used to study the tunneling 
dynamics of ultracold bosons in a double well potential in [3, Q and could be used also for Fermi gases. The energy 
barrier and the distance between the wells are controlled by acting on the parameters of the optical lattice |11| , while 
the number of doublets under barrier depends also on the trapping frequencies of the harmonic confinement. In the 
following we discuss how the fermions fill the doublet structure, keeping into account the transverse eigenstates. We 
will also study two relevant situations: the ID limit with frozen transverse degrees of freedom and the limit (which 
we refer as 3D) in which the transverse frequency is so small that only few doublets of the double well potential are 
occupied. 

The trapping potential reads 

V(x.) = ~mul(x 2 + y 2 ) + V DW (z), (1) 

where the double well potential along the z-axis, Vbw( z )j has the form 

Vdw(z) = -zmuj\z 2 + Vo cos 2 (kz). (2) 

In ©, the ID optical lattice is created by two counterpropagating laser beams: it is k = 27r/A where A = 
Xinx p. r sin (0/2), Xi aser being the wavelength of the lasers and 9 the angle between the counterpropagating laser beams 
[ll| . The spacing in the lattice is A/2, and Vq is proportional to the power of the laser; moreover, oj z and u>± are 
respectively the axial and transverse frequencies of the harmonic confinement. Typical experimental numbers are 
A ~ 1 — lOfim and uj z /2tt ~ 10 — 100Hz: for A ~ 10/im, the energy barrier (needed to have tunneling dynamics) is 

Vo/h > 500-ffz 7]. In the following we set Vb = sEr, where En — 4p- is the recoil energy, and we also introduce the 
ratio 

A = ^ (3) 
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FIG. 1: Solid line: plot of the potential Vdw(z). We used the values cj z = 2-7T-78 Hz, A = 10.4/im and Vb/ft = 13E_r//i = 3.7kHz 
corresponding to s = 13 and A = 0.17. Dashed lines: harmonic potential around the minima of the original potential with 
effective frequency ftw — 7.2Er. 



where auo = ^h/muj z is the harmonic oscillator length in the z-direction. 

For suitable ranges of the parameters, the potential ([1} is to a very good approximation a double well potential, 
as it is shown in FigJT] where we fixed the value of the frequency of the harmonic trapping along the z-direction. The 
potential Vdw(z), rescaled in terms of Er, reads: 



Vpwjz) 
Er 



= s cos 2 (kz) + 



(2ttA)' 



-(kz) 



(4) 



The first two minima of the potential which identify the position of the double well minima are given by the solution 
of the transcendental equation 



sin (2kz) 

2kz 
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(5) 



One sees that for large values of the denominator in the r.h.s. of Eq.Q - i.e. for large energy barriers, s>l- the first 
two minima are very close to kz = ±5. Expanding the potential around these minima we get the effective trapping 
frequency 



2y/sE t 



(6) 



the energy difference between (neighbour in energy) doublets is then ~ Huj. This frequency then allows to estimate the 
mean energy and the number of doublets under the potential barrier Vq ■ In general the splitting of doublets increases 
when the energy of the doublets increases. A first estimation of the splitting of each doublet can be performed by 
means of a semiclassical computation |37j |: in Appendix we report semiclassical estimations of the energy and the 
splitting of the doublets. 

In Fig[2] we plot the potential and the position of the doublets for specific values of the parameters. In Fig(3] we 
plot the number of doublets for the potential of Eq.(H} for different values of Vq (in units of Hu) z ) and A = 2ir/k 
obtained from the solution of the time-independent Schrodinger equation. We observe that for a fixed ratio a^ol 'A 
there is a finite range of values of Vq' for Vb smaller than a critical value y^ mm > the central barrier is not high enough 
to contain any doublet. For large values of Vq, the nearest energy minima (located, for s ^> 1, close to kz — ±4^) 
have an energy smaller than the potential energy at the top of the barrier, which is sEr, and then the double well 



structure disappears. From the condition Vdw(z 
given by 



±|?) > Vdw(0) one gets the maximum value y^ max > ^ which is 



V, 



(max) 



9 

64tt 2 A 4 



En 



32 



;?A 2 . 



(7) 




FIG. 2: Plot of the potential Vdw(z). Here we use the values uj z — 2% ■ 78 Hz, A = 15/itm and Vo/h — 35E_r//i = 4.8kHz 
corresponding to s = 35 and A = 0.12. The double well has three doublets under the central barrier; the upper solid line 
correspond to the symmetric energetic level of the third doublet; the dashed one to the antisymmetric level. For the first two 
doublets (solid lines) it is not possible to distinguish among symmetric and antisymmetric levels because are too close in energy. 




FIG. 3: Number of doublets as a function of the strength Vb and the wavelength of the periodic potential expressed, respectively, 
in terms of hv z and auo- 



A. Number of particles and dimensionality 



Our goal is to study the dynamics with an initial condition characterized by an imbalance in the number of atoms 
between the two wells. To this aim, we have to determine first how the fermions fill the energy levels. Since typically 
Q ^> u> z , then the system is essentially ID when 



Nw<^uj_ l , 



(8) 



where N is the total number of particles. In this case fermions begin to occupy the levels of the double well potential 
along z; the number of fermions is given by two times the number of doublets. 
The opposite 3D limit is obtained when 



Nuj± < LU. 



(9) 



In this limit, all the particles stay in the first doublet, occupying the transverse eigenstates associated to the doublet 
eigenstates. More precisely, to accommodate all the particles in the first doublet, two conditions must be satisfied 
at the same time: the energy difference between the first and the second doublet states (~ luij) must be much larger 
than the energy of the transverse modes (ex Huj±) and the transverse trapping frequency uj± must be much smaller 
than uj z . Indeed, if oj± ~ ui z , then the two conditions u = ky/Vo/m 3> ijJ± ~ u> z and Vo < y^ max > — (9/32)mw 2 A 2 
(needed to have a double well structure) cannot be at the same satisfied. Finally, we notice that if Nuj± > uj then 
higher doublets can be occupied and situations intermediate between the ID and 3D limits can be explored. 
Before studying the 3D double well potential ([T]), let us consider a 2D potential of the form 

V 2 d(x,z) = -mujx 2 + -mujlz 2 + V cos 2 (kz). (10) 

The number of particles which can be placed on a side of the double well in the first doublet is then given by 

N 2D = ^ (11) 

nuju 



where 



a Ei = (l±4hM±^l (12) 



is the energy difference between the first and the second doublet: e 1 \ is the energy of a symmetric (S) or antisymmetric 
(A) state of the first (1) or second (2) doublet. For very large barriers, then AEi ~ Mi [37(. If the energy of a state 
in the second doublet is larger than the sEr than the number of particles that it is possible to store up to the barrier 
height is given by: 

sE R 



N 2D = I= « ^_ (13) 

own 

In three dimensions the number of particles which can be stored is even larger thanks to the degeneracy of the 

transverse degrees of freedom. The energy of a state in a side of the barrier is given by hu!±(n x + n y + 1) -\ — " 2 " , 
where n x ,n y = 0, 1, • • • denote the transverse quantum numbers and n — 1, 2, • • • denote the doublet. If we denote 
by h = n x + n y then the order of degeneracy of the transverse level is given by fi + 1 whereas the total number of 

levels up to h is given by — — -^ — - ~ \. Then the total number of particles which can be stored on a side of the 
double well in the first doublet is given by: 

an 1 f*Ei 



NW = 2 ( >srJ ■ (14) 

When the energy of the second doublet is larger than sEr then as for the 2D case, AEi must be replaced by sEr — 

e s +e A (e S +e A )~(e s +e A ) 

1 2 1 . These results can be extended to higher doublets by simply substituting AE\ by AE n — " + — ^^- — — — —. 

III. TUNNELING DYNAMICS IN A DOUBLE WELL POTENTIAL: FORMALISM 

In this Section we review the formalism to describe the dynamics of a system of polarized fermions in a 3D double 
well potential of the form Eq.(fT]). The ID limit, studied in [34|, is retrieved and briefly discussed; we refer also to [34[ 
for a study of the effects of the boson-fermion interactions when a localized Bose-Einstein condensates is placed in 
the wells. 

The Hamiltonian for a system of polarized fermions confined in an external potential V(x) reads 



H = dx ^(x) 



-^V 2 + y(x) 
2m 



V(x). (15) 



The eigenvectors of the single-particle Hamiltonian 

h 2 



-2^ V2 + y W 



<PnAx)Vn y {y)4> 1 {z) = e nxnyl ip nx (x)ip ny {y)^(z) (16) 



are defined by three quantum numbers: n x , n y = 0, 1, • • ■ denote the transverse quantum numbers associated to the 
transverse harmonic oscillator, while 7 denotes the eigenstates of the double well potential Vr>w{ z )i as obtained from 



2 n2 



h 2 & 

2m dz 2 



Vdw{z) 



4> 7 (z) = e 7 (/> 7 (z). 



(17) 



The quantum number 7 is defined by the pair n, a, where n — 1, 2, • • • denotes the doublet, and a — S,A denotes 
the symmetric or antisymmetric state (in the following we denote e 7 by e"). The eigenenergies of the single-particle 
Hamiltonian are given by e„ xrl 7 = hujj_(n x + n y + 1) + e 7 . 

The fermionic field can be written as ^(x) = J~] ip na: (x)ip ny (y)(j) 7 (z)c na . nyl . Since we want to study the 

dynamics of particles in the double well potential we decompose each pair of eigenfunctions corresponding to the 
same doublet into its right and left component: 



^( X ) = Yl VnA^nyiy) 4>n{z)Cn x n y n + (t>n{ z ) c 






(18) 



where the (Wannier) wavefunctions 4> R ' L centered in the wells R, L are given as usual by 

<f>n,s{z) ~ <t>n,A{z) 



k R( ^ _ <P n,s(z) + (f>nA Z ) ,L( 

' " 71 ' K{ *> ~ 7T 



CW 



We will as well use the notations 4> nic n y n,a{x-) — { Pn x { x ) i Pn y {y)4'n,a{z) for the eigenfunctions of the single-particle 
Hamiltonian and ip^'n ra ( x ) — L Pn I {x)ip n (y)<fi R ' L (z) for the 3D Wannier wavefunctions. The fermionic operators 
c R ^ „ arc defined as: 



"n x n y n,S < Cn x n y n.,A j^ 



S - Cn 



-n x n y n^ ^n^n^n 



V2 



"> ^n x n v n 
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The Hamiltonian (TT5t then reads 



e s + ( A 



H = ^ <| I hw\_ (n x +n y + 1) + , } 
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(19) 



(20) 



9 \ n x n y n L 'n x n y n ~ ^ J rixn 1J n L 'n x n y n 



The density operator of particles along the direction of the double well is given by: 

n{z) = J dxdy V> f (x)^(x) = ^ ^ 7 (z)0y (z)c^„ H7 c„ x „ y 

n. ™ n. „ . -v 0/ ' 



(21) 



n.n, 7l7 ' 



When we average the density operator over a generic state with a definite number of particles on the left and right 
part of the barrier the off-diagonal elements of the density operator give a null contribution. Then, in order to describe 
the evolution of the system, we write the equation for the dynamics of the diagonal elements only: 
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(22) 



and similarly for c R '^ n c^' R n- Finally one easily obtains: 



by It, I b x I by 
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(23) 



We introduce the quantity 



A7V (t) = /c flt c R \ 



( c^„ n c L n „ „ 



which is the particle difference in the doublet n x ,n y ,n between the right and the left sides of the double well for a 
generic initial state -00 (given an operator 0, we use the notation (O) = (ijjq \ 0(t) \ ipo))- Next we define the total 
population imbalance AN n (t) for the n-th doublet: 

AN n (t) = J2 ((<£[»,»<&»,«) ~ (<#„,»<£„„,„)) • (24) 



The total fractional population imbalance z(t) is then defined as 

E„ajv„(<) 



z{t) = ^" N ny ' - (25) 



From Eqs. ([2"3")) . AN nxn n (t) evolves according to: 

^AN nscnyn {t) = (£j^j AN nocnyn (t). (26) 

This equation has the straightforward solution: 

AN nasnyn {t) = AN nxHyn {0) cos (2Q n t + </>„,„„„) , (27) 

e A -e s 

where fl n — n 2h " is the Rabi frequency for the oscillation of a particle into the same doublet and ip nn .n y n depend in 
general on the initial conditions. 

For a single particle in the doublet n x , n y , n, the initial state ^o( x ) can be in general be written in the Bloch sphere 
as 

V>o(x) = cos f " x " y " J </>„ x n y n,s(x) + e^ 71 -"*™ sin ( "*""" J (^^^(x), (28) 

where 6*„ xny „ and <fn x n y n are the coordinates on the Bloch sphere: when the particle is initially in the state <^ „ n (x) 
( < /'™ x n ! ,n(x)), then #„ x „ H „ = 7r/2 and ip nas n v n = (<y5n,«„n = t). The energy of such a state is equal to: 

{MH\i> ) = 4 cos 2 (^f^j + ei sin 2 (^f^j . (29) 

It is then easy to see that the population imbalance for this state evolves according to (|2"T|) where 

AN n:cnyn (0) = sme nxnyn . (30) 

Let consider now many fermions in the double well potential. For a generic initial state with the fermions in the same 
n-th doublet with different initial phases y na .n n, then oscillations in the population imbalance are rapidly washed 
out: indeed, since the fermions are polarized, the dynamics in different doublets are independent and dephasing of the 
different particle differences AN n!cnyn (t) in general occurs. However the experimental situation we have in mind is a 
system of two identical wells initially practically decoupled (i.e., with a very high energy barrier between them) with 
a definite number of particles (say Nl and Nr, with Nl < Nr) in each side in their ground state. Then, at the initial 
time t = the barrier along the z direction is lowered and the system is left to evolve freely with the hamiltonian (|15l) . 
Alternatively, one could start from the ground state of a symmetric 3D double well and tilt for some time the double 
well, so that the energy minima of the double well potential are different and particles flow toward the energetically 
favoured well: removing abruptly the tilted potential one has an initial state with a different number of particles in 
the two wells. 

For sufficiently high barriers the first Nl particles will be frozen in the two wells, while the others will start to 
tunnel. Take as initial state for one of such particles in the right side ijj nx n n z where n x and n y are the quantum 
numbers along the transverse directions which remain unperturbed and n z labels the states along the longitudinal 
z-axis before the barrier gets lowered; then the initial state t/>o( x ) = Vv^n n z (x) can be well approximated by 

V>o(x) « <^„ B „(x). (31) 

Therefore we decide to study the evolution of the system of an arbitrary number of fermions under the conditions 
that the initial state of each particle is of the type ip nx (x)ip ny (y) J- n to have a reasonable description of the 

dynamics. 
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FIG. 4: In a), b) c) we plot z(t) for a ID double well ([2]): in a) there are 4 fermions initially on the right side and 3 on the 
left side, in b) 5 on the right and 3 on the left and in c) 6 on the right and 1 on the left. The parameters of the potential are: 
Vo = 20huz and A = 15aho- In d) it is plotted z(t) for the square double well potential (|34[) with 100 fermions initially on the 
right and 60 on the left. Parameters: a = 1, b = 5 and Vo = 3200 (with h = m — L — 1); notice that time is in logarithmic 
scale. 



If, for instance, we suppose all the particles in the n-th doublet initially on the right well the phases </?„ x „ „ can be 
taken equal to zero and the number of particles in each side of the barrier read: 



* = J2 (cn!n y nCtn y n) = AJV„(0) CO S 2 (O n t) 



(32) 
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A. Dynamics in the ID limit 

In the ID limit the transverse degrees of freedom are frozen and fermions belonging to different doublets oscillates 
with different frequencies: this case been studied in detail in J34j| . In Fig0]we plot the evolution of the total fractional 
population imbalance z(t) for few particles (N = 7 and N = 8). In the last plot we draw the same quantity for 
N = 160 in a square double well potential (we made this choice since for the potential like (|4]) one cannot have 
hundreds of doublets below the barrier for realistic values of the parameters). We parametrized the square double 
well potential with infinite walls as 



V(x) 



2mL 2 





V \x\ < a L 
a L < \x\ < b L 

d \x\ > b L 



(34) 



where 1L is the width of the square double well and a, b dimensionless parameters. One sees that already for hundreds 
of particles the current oscillations are practically washed out for an initial imbalance not too small. 

As we can infer from these examples, oscillations follow an aperiodic motion due to the superposition of many 
different incommensurate oscillations; moreover, since Rabi frequencies for distant levels can be very different, the 
time scale required for a complete oscillation of the population imbalance z(t) can be much larger than the experimental 
observation times. 




FIG. 5: Plot of the maximum number of particles which can occupy the first doublet as a function of the strength of the optical 
lattice potential Eq.© with u± = 10Hz, v z — 300Hz and A = 0.25. For these values of parameters, there is only a doublet 
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<V < V^ max) with V y 
' 3 (circle) and Vo/frw z 



' (mm) /^ = 
Vo/Kuz = 7.5, N = 5.1 ■ 10 a (circle) and V /hu z = 7.5, N = 
occupied and the second one is almost completely occupied. 



2.1 and Vq /fajj z — 4.5. The circle and the star correspond to 
= 5.1 ■ 10 3 (triangle): in both cases the first doublet is completely 



IV. DYNAMICS IN A THREE-DIMENSIONAL DOUBLE WELL POTENTIAL 



In this Section we consider the dynamics of polarized fermions in the 3D confining potential ([T]) when the ID validity 
condition is violated, as it happens for realistic transverse confining potentials. 

As discussed in Section II, the number of particles which can be stored in one doublet can be very much increased 
by taking a weak confinement in the transverse direction thanks to the possibility to fill the transverse states before 



reaching the next doublet. In FigJS]we consider the range of parameters Vq < y Q < yv max ) gucn that the double 

well potential Q has only a doublet under barrier and we plot in this range the maximum number N max of particles 
which can be stored in that (first) doublet: in other words, for a number of atoms N > N max they start to occupy 
levels above the top of the barrier. One sees that N max ~ 10 3 for reasonable experimental parameters of the potential. 
Since in this case there is only one doublet with energy below Vb then the number of particles scales as Vq (the energy 

of the particle in the double well direction, 1 1 , grows as VVo ) ■ 



The number of particles can be further increased by decreasing the ratio A = ^- , as shown in FigO In this figure 
we consider values of the parameters such that the double well potential (j4]) can have one, two or three doublets: we 
plot the maximum number of particles N n which can be stored in the doublet n — 1, 2 or 3. The maximum number 
of particles in the first level can arrive to 2V max « 8 • 10 3 . 

The dynamics for the fractional population imbalance z(t) is illustrated in FigjTJ here for simplicity we show the 
results for the evolution of particles initially distributed in the the right side of the barrier. In the left part of the 
figure we considered N = 5.1 • 10 3 atoms distributed among the first two doublets filling completely the first one and 
almost totally also the second at Vb = 7.5Huj z (corresponding to the circle in FigJS]): we see that the oscillation of z(i) 
contains only two Rabi frequencies. The (normalized) particle density at different times is plotted in FigfS] In the 
right part of Fig[7]we considered N = 9.7 • 10 3 atoms distributed among the first (completely filled) and the second 
doublet at Vb = 10huj z (corresponding to the triangle in FigJS]), and the particle density at different times is plotted in 
FiglHl In this case, the profile of the densities has more pronounced secondary peaks due to the contribution coming 
from particles in the second doublet. 

Finally we observe that a way to observe the relative number of particles in each doublet at the initial time t = 
is through the analysis of the Fourier transform of the fractional population imbalance. Since our considerations are 
done at T = the motion of the particles in the system is undamped, so the Fourier transform of z(t) is simply 
given by a sum of (^-functions, each one localized at the characteristic Rabi-frequency of each doublet: for the cases 



10 




5 10 

VJhv 

z 



FIG. 6: Plot of the maximum number of particles Ni, N% and N3 which can be stored respectively in the first, second and 
third doublet as a function of the strength of the optical lattice potential in Eq.@ with u± = 10Hz v z — 300// z and A — 0.15. 
For these values of parameters up to three doublets have energy lower than Vo . 




FIG. 7: Above: (full line) evolution of z(t) for the points of Fig[6] the circle (left) and the up triangle (right); (dashed line) 
evolution for a single particle in the lowest doublet with equal initial conditions. Below: magnitude of the modulus c(w) of the 
Fourier transform of z(t) for the two plots above. 



considered in Figs|BllHl where only two doublets are occupied, one has 

r- +00 



z(u) 



e wt z{t)dt 



N 9 



2,v 1 ^ [s(u - no + s(lu + no] + ^ [s(u - n 2 ) + s{lu + n 2 )] 



(35) 
(36) 



where N\ and N2 are the number of particles respectively in the first and second doublet. In the inset of FigjT] we 
plot the coefficients of the Fourier transform for the two cases described above. 
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FIG. 8: From above left: normalized particle density at different times for the circle of Fig[6] We take 5.1 ■ 10 3 particles 
distributed among the first two doublets, a), b), c), d): ui z t — 0, 100, 225, 400. 




FIG. 9: From above left: normalized particle density at different times for the triangle of FigEJ We take N = 9.7 ■ 10 3 particles 
distributed among the first two doublets, a), b), c), d): u>zt — 0, 250, 750, 1000. 



CONCLUSION 



We studied the tunneling dynamics of a spin polarized Fermi gas in a three-dimensional double well potential at 
zero temperature. We focused in particular on the time dynamics starting from an initial state in which there is an 
imbalance between the number of particles in the two wells. Although fermions in different doublets of the double 
well tunnel with different frequencies, we point out that (incoherent) oscillations of a large number of particles can 
arise, as a consequence of the presence of transverse degrees of freedom. 

Estimates of the doublet structure and of the occupation of transverse eigenstates for a realistic experimental 
setup are provided. In the ID limit the current oscillations are washed out, as a result of the dephasing, but for 
not too large confining transverse frequencies or not too small number of particles the fermions can occupy only the 
first doublets (using the transverse eigenstates), resulting in oscillations of a large number of fermions. We stress 
that these are incoherent single-particle oscillations. We can conclude that if the ID condition is violated, then the 
simple observation of a sinusoidal current cannot in general discriminate between incoherent and coherent ( Josephson) 
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Doublet 


A-Efixact 


Ai?Scmi 




1 
2 
3 

4 


5.44 • 10" b 

3.34 ■ 10" 3 

8.75 ■ 10" 2 

1.12 


4.95 • 10 -5 
3.70 ■ 10~ 3 

0.142 

2.80 


-0.09 
0.11 
0.62 
1.50 



TABLE I: Semiclassical estimation of the energy using the value given by (JB]) for the energy E. Energies are in units of En. 
We use here a = 76 and A — 0.1 corresponding to four doublets. 



Doublet 


Ai^Exact 


A^Semi 


EScmi--E E xact 


1 
2 
3 

4 


5.44 • 10 _b 

3.34 • 10 -3 

8.75 ■ 10" 2 

1.12 


4.93 • 10 -5 

3.19 -10 -3 

8.85 ■ 10 -2 

1.14 


-0.094 

-0.045 

0.011 

0.018 



TABLE II: Semiclassical estimation of the energy using the numerical value for the average energy E. Energies are in units of 
E R . We use s = 76 and A = 0.1 as in TableU 

tunneling dynamics at zero temperature. 
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Appendix A: Semiclassical estimation of the splittings 



In this Appendix we report a semiclassical estimation of the energy splitting among the levels of the double well 
potential ^. 

The semiclassical formula for the splitting of the n-th energy level is given by |37| 






hid 



± P -\r-o,^ \p\ 



(Al) 



where a,~a denote the turning points at the central barrier of the potential (classical motion is inhibited there), 
\p\ = y // 2m(V(z) — E) and w c i is the frequency of the classical motion between the turning points where classical 
motion is allowed. 

In the Tab.(|lJ we compute the semiclassical splitting energy using the estimation for the average energy E of the 
doublet: 



E = lil' \n+ - 



(A2) 



based on the effective frequency (|6]). 

In the Tab.(jn| we compute the splitting energy using the average value of the energy of each doublet obtained from 
the numerical solution of the Schrodinger equation for the potential ([2]). 

As we might expect the agreement of the semiclassical estimation with the exact values is much better in the second 
case; we also note that in the first one there is a good agreement (the error is around 10%) for the first two doublets 
due to the validity of using the effective frequency (j5]) for the lowest energy levels. 
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